A local Monte Carlo implementation of the non-abelian 
Landau-Pomerantschuk-Migdal effect 
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The non-abelian Landau-Pomeranschuk-Migdal (LPM) effect arises from the quantum interfer- 
ence between spatiaUy separated, inelastic radiation processes in matter. A consistent probabilistic 
implementation of this LPM effect is a prerequisite for extending the use of Monte Carlo (MC) 
event generators to the simulation of jet-like multi-particle final states in nuclear collisions. Here, 
we propose a local MC algorithm, which is based solely on relating the LPM effect to the prob- 
abilistic concept of formation time for virtual quanta. We demonstrate that this implementation 
of formation time physics alone accounts probabilistically for all analytically known features of the 
non-abelian LPM-effect, including the characteristic L^-dependence of average parton energy loss 
and the characteristic modification of the gluon energy distribution. Additional kinematic 

constraints are found to modify these L^- and tj-dependencies characteristically in accordance with 
analytical estimates. 



In heavy ion coUisions, highly energetic partons prop- 
agate through the produced dense QCD matter before 
reaching the QCD vacuum. The fragmentation of these 
partons is modified significantly by the medium. This is 
seen in the strong suppression of essentially all hadronic 
spectra at high transverse momentum, discovered at 
RHIC T|, and is consistent with calculations of medium- 
induced parton energy loss 0, H, 0, @, S Q- Heavy 
ion collisions at the LHC will open qualitatively novel 
opportunities for studying this jet quenching effect Q, 
since jet-like multi-hadron final states will become exper- 
imentally accessible in a wide kinematical regime, which 
is well-separated from the soft high-multiplicity features 
of nuclear collisions. To exploit these novel opportuni- 
ties fully, a reliable theoretical framework for the de- 
scription of multi-particle final states and their medium- 
modification is needed. In the absence of a medium, the 
standard tool is the MC simulation of the final state par- 
ton shower underlying jet fragmentation. Its dominant 
medium-modification is expected to be determined by 
the non-abelian LPM effect. The question arises, how 
this medium-induced quantum interference effect can be 
implemented in a probabilistic MC simulation. Here, we 
propose a MC algorithm to this end. 

In the vacuum, a virtual parton p — q, g of energy 
E degrades its virtuality mainly by radiating a gluon 
of energy uj and transverse momentum k±. The gluon 
energy distribution is £^ ^ Pg/piz) where 

Pg/p{z) ^ ^ + 0(2°), z — uj/E denotes the standard 
p p + g splitting function in the vacuum. The 
consequences of embedding the parton splitting process 
p —^p + g into a medium have been studied in de- 
tail [a H, S, 0, 0, Q- For relativistic energy, the lead- 
ing medium-modification can be written in terms of an 



Pomerantschuk-Migdal (LPM) effect. This contribution 
persists up to an energy scale w < Wc = ^9^^, defined in 
terms of the in-medium path length L and the quenching 
parameter q @, 0, 0] . 

There have been several MC approaches to the prob- 
lem of radiative parton energy loss. These implement the 
LPM effect for instance by some form of parametrized 
suppression of the projectile energy or a parametrized 
modification of the parton splitting 0, [ll|. Destruc- 
tive interference in inelastic processes has also been mod- 
eled by assigning to newly produced partons a mean 
free path, during which no further interactions can oc- 
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14| . Moreover, medium-induced parton energy loss 
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which shows the char- 



acteristic Xj^fuj- dependence of the non-abelian Landau- 



has been recognized to change the virtuality evolution 
of the parton shower (see e.g. [1, [3]) and there is a MC 
model implementing a medium-modified virtuality evolu- 
tion The MC algorithm proposed here differs from 
these approaches in that the LPM effect does not emerge 
as a consequence of a parametrized modification of the 
parton cascade. Rather, it is the dynamical result of 
a local, energy-momentum conserving evolution, which 
is consistent with analytically known limiting cases of 

To understand how LPM quantum interference can 
be built up in a local probabilistic dynamics involving 
multiple scattering, we recall first analytical results in 
known limiting cases We start with the case of a 
gluon radiation process involving two transverse momen- 
tum transfers q\^,\ and q^.2 within path length L. The 
resulting distribution of radiated gluons is known ana- 
lytically ioY E ^ u: ^ \k\_\. It interpolates between 
two limiting cases. In the totally incoherent production 
limit, the gluon is produced in the first scattering and its 
transverse momentum is shifted probabilistically by g_L,2 
in the second scattering. In the opposite limit, the two 
scattering centers are not resolved and the radiated gluon 
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spectrum corresponds to that of a single (effective) scat- 
tering center, in which the momentum {q±^i + q±,2) was 
exchanged with the target The transition between 
these limits occurs at a scale L ~ l/Qi = 
where is the transverse momentum of the radiated 
gluon after the first momentum transfer. This scale 
enters the analytical result via the interference factor 
2 (1 — cosQi L) I {Qi L)^, which vanishes in the incoher- 
ent limit L 3> l/Qi and which approximates unity in 
the totally coherent hmit L < l/Qi Q- The term l/Qi 
can be interpreted as the gluon formation time after the 
first scattering. If the second scattering occurs within the 
gluon formation time, then it acts coherently, otherwise 
incoherently. By inspecting analytical results for gluon 
production processes involving more than two momen- 
tum transfers [3] , one finds that this role of the inverse 
transverse gluon energy as formation time is generic [l3| : 
the next (n -I- l)-th scattering center contributes totally 
coherently (incoherently) to a gluon production process, 
if the gluon formation time 
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after the n-th scattering is much larger (much shorter) 
than the distance d to the {n + l)-th scattering center. 

These limiting cases of totally coherent and incoherent 
gluon production are implemented easily in a probabilis- 
tic MC algorithm. Consider a target scattering center 
Qt participating in an inelastic process. In the incoher- 
ent case, the gluon will be produced on this scattering 
center with inelastic cross section (jiQT^iQTg ^ Subse- 
quent scatterings can change with elastic cross section 
(jgQr^gQT ^j^g transverse momentum of the fully formed 
gluon to kx ^ kx + qA_.i, but they are not necessary 
to complete the gluon radiation process. In the opposite 
limit of totally coherent gluon production, all n scattering 
centers within the path-length L transfer their individual 
momenta q±^i coherently to the projectile, as if the gluon 
were radiated in a single inelastic scattering process with 
transverse momentum transfer 



q±,tot 
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Although this limit is dominated by quantum interfer- 
ence, it can be realized exactly in a MC algorithm by 
selecting probabilistically those target scattering centers, 
which transfer a non-vanishing transverse momentum 
q±^i from the medium to the projectile. The distribu- 
tion of the q±^i is determined by the elastic cross section. 
The gluon is then generated according to the cross section 
g-qQr^qQTQ with the constraint that the momentum 
transfer from the target is q±^tot- This procedure matches 
exactly the result of analytical calculations [1] . 

In between the totally coherent and incoherent limits, 
we propose that the LPM effect can be implemented in 



a probabilistic MC algorithm by requiring that the mo- 
mentum transfer from different scattering centers to the 
partonic projectile acts totally coherently for gluon pro- 
duction, if it occurs within the formation time tf, and 
that it acts incoherently, if it occurs after tf. This is 
implemented as follows: Once a scattering center is se- 
lected as the source of an inelastic gluon production pro- 
cess, the MC determines probabilistically the distance d 
to the next scattering center with which the projectile 
(quark and gluon) interacts. If d > tf, then the gluon 
is regarded as formed, and it is treated in subsequent 
steps of the dynamical evolution as an independent pro- 
jectile. If d < tf, then the momentum transfer of the 
subsequent scattering center is added coherently to the 
one in the gluon production process. This changes the 
gluon transverse momentum, and thus the gluon forma- 
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Then, if the next interaction 



lies outside the length t'j, the gluon is regarded as formed. 
Otherwise, one repeats the above procedure, until the 
gluon is formed. Another inelastic interaction is only 
generated once the gluon is fully formed. This defines a 
local, probabilistic time evolution, which reproduces by 
construction the analytically known totally coherent and 
incoherent limiting cases. 

We now demonstrate in a numerical model study that 
such a MC algorithm reproduces the main character- 
istic features of the LPM effect. To be specific, we 
model the target as an ensemble of colored scatter- 
ing centers Qt with density n. The scattering centers 
present elastic and inelastic cross sections (jiQi'^'iQT and 
^qQT~*qQTg projectile parton. Here, we use the 

ansatz A-ai^T^qQTg = g2 (j^, 1 fjqQr^qQr ^ which was 
used in the calculations of [1, [1^ 0, d, H, 01 • The results 
shown below were obtained for total elastic and total 
inelastic cross sections fixed such that n cr'^'^^^'^'^^s = 
Ifm""'^ and ria^'^'^^^''^'^ = lfm~^. This correspond to 
q ~ IGeV^/fm, where q characterizes the strength of 
elastic scatterings 
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dfjqQr^qQi 
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We checked that consistent with analytical results 0], 
the q- and i-dependence of the MC simulated spectrum 
'^^'^ accommodated in a function of lu/lOc (data 
not shown). As a consequence, our numerical results do 
not rely on a particular parameter choice for q and we 
plot them in suitably rescaled diniensionless quantities. 

We have calculated the average energy loss AE' for a 
projectile quark of initial energy Eq = 100 GeV travers- 
ing an in-medium path length L. Results for the full 
MC algorithm, including LPM-effect and exact energy- 
momentum conservation at each vertex are shown as 
' full case r in Fig.[TJ All simulations are done for an elas- 



tic cross section 



W+q^ with /i = 1 GeV 



We recall that the BDMPS multiple soft scattering ap- 
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FIG. 1: The medium-induced energy loss AiJ of a quark of 
initial energy Eq — 100 GeV as a function of the in-medium 
path length L. Results for the full MC algorithm with LPM 
effect and energy conservation (full case 1) are compared to: 
i) LPM effect without exact energy conservation (case 2), ii) 
incoherent limit without energy conservation (case 3), iii) in- 
coherent limit with energy conservation (case 4), and iv) full 
case 1 with constraint that momentum transfer per scattering 
center is limited to q± < 2^ (soft case 1). The bottom plot 
zooms into the small L-region. 



proximation neglects the high-g^ tail of this distribu- 
tion To mimic this approximation, we consider also 
a ' soft case 1 ' , in which the transverse momentum range 
is (7_L G [0,2^] rather than e [0,Ci;]. 

For small L, results for the soft case 1 show an exact 
A.E cx L^-behavior (see Fig.[T)3), and the gluon spectrum 
shows the characteristic ~ 1/ y/ui. [We find dl/dui cx 
l/w^-^^ for u! < 30 GeV, data not shown.] This demon- 
strates that the MC algorithm reproduces the character- 
istic BDMPS results if the soft BDMPS approximation 
is implemented. If one removes this soft approximation, 
q± £ [0,w], one finds AE oc L^'^, which lies in between 
the quadratic dependence of coherent and the linear de- 
pendence of incoherent radiation. This is so since large 
momentum transfers > 2fi can shorten the formation 
time abruptly, thus introducing for small L an incoherent 
component of gluon radiation. For the same reason, we 
find for the full case 1 and w < 30 GeV a best fit to the 



spectrum dl/duj cx l/oj^ , where /? = 1.38 lies slightly fur- 
ther away from the BDMPS-limit {(3 = 3/2), and closer 
to the incoherent limit {(3 — 1), see Fig. [2l 

For sufficiently large L > Lc — \J^^-^ a complete 
break-down of coherence is expected to occur, since the 
kinematic constraint uj < Wmax implies a bound on the 
average formation time (t/) < 2ujjnax/ {k'^)- Within seg- 
ments of length L < Lc, gluons can be emitted in a 
completely coherent process. But different segments of 
length Lc add incoherently, and hence AE turns from a 
quadratic to a linear dependence around L ^ Lc, consis- 
tent with As AE grows with L, energy conservation 
becomes more important and AE approaches Eq for very 
large L, see Fig. [T] 

In case 2 . we study the high-energy approximation that 
the projectile energy is so large that its reduction due 
to gluon emission can be neglected. For small L, when 
energy conservation is unimportant, results on average 
energy loss and the gluon distribution agree with case 
1, see Fig. [T}) and [5^. For large L, however, since the 
projectile energy is kept fixed by assumption, all seg- 
ments of length Lc contribute the same, and AE cx L as 
L ^ Lc (see Fig. [T^). Moreover, as for the incoherent 
limit, dl / dhj ~ l/oj for very large path length L/Lc 3> 1 
(data not shown), but the gluon yield remains strongly 
suppressed in comparison to the incoherent limit, since 
scattering centers within a distance Lc act coherently. 

We finally compare these results to the incoherent lim- 
its case 3 and case 4 . obtained by setting tf = in all 
intermediate steps of the MC algorithm. In case 3, the 
energy of the projectile is kept fixed as in case 2, i.e. ex- 
act energy conservation is abandoned, but one requires 
Lo < Winax = Eq. The average energy loss is then linear 
for arbitrary path length L (see Fig.[T]). For our model of 
the target, dl/dui ^ I/lo (see Fig. [2]), and the slope of the 
linear dependence of AE cx L depends logarithmically on 
the UV-cut-off Wmax- In case 4, we repeat the calculation 
with exact energy conservation. For small L, this con- 
straint is negligible, and results agree with those of case 
3. For large L, however, energy loss is dominated by en- 
ergy conservation and AE approaches Eq. Also, energy 
conservation implies that at large L, only soft gluons can 
be emitted, since the projectile has lost already most of 
its energy. As a consequence, the spectrum at large L 
tends to be much softer than dl/dcu cx l/us. [The best 
fit in Fig. [^b is dl/doj cx l/w^'^, remarkably close to the 
LPM-spectrum.] 

To sum up: In vacuum, the dominant quantum in- 
terference between subsequent parton branchings can be 
taken into account probabilistically by angular ordering. 
Here, we have shown that the corresponding dominant 
medium-induced quantum interference can be taken into 
account probabilistically by a formation time constraint. 
In particular, this reproduces by construction the ana- 
lytically known limiting cases of totally coherent and in- 
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FIG. 2: The medium-induced gluon energy distribution radi- 
ated from a quark of energy Eq = 100 GeV for the four cases 
of Fig. [Hand path length (a) L/Lc ~ 0.5, (b) L/L, ^ 2.5. 



coherent gluon production, and it reproduces the char- 
acteristic AE cx L^- and ujdi/duj cx l/-y/aJ-dependencies 
of the non-abehan LPM effect, which were calculated in 
the soft scattering approximation i (soft case 1), and 
it implements the known scale of breakdown of coher- 
ence effects ^] at Lc- By comparing simulations with 
medium-induced quantum interference to the incoherent 
limits (cases 3 and 4), our study also further quantifies 
the known importance of the LPM-effect in the calcula- 
tion of radiative parton energy loss. 

Reproduction of known analytical limits is an impor- 
tant test for a MC algorithm of the LPM-effect. Finally, 
however, a major use of MC algorithms lies in exploring 
physics effects, which are difficult to treat in analytical 
formulations. The proposed MC algorithm is suited to 
extend analytical treatments in several directions: First, 
the analytical results of [3, [3, 0, 0, 0, 0] rely on approx- 
imations, in which energy is not conserved exactly, and 
the energy of the incoming projectile remains unchanged 
in subsequent gluon emissions. In the present study, this 
corresponds to case 2, and it is seen to differ significantly 
from a formulation with exact energy-momentum conser- 
vation (case 1) outside a region of relatively small path 
length and very high projectile energy. Second, the ana- 



lytical calculations of 0, 0, 0, H, S 01 do not consider the 
possibility that emitted gluons undergo further medium- 
induced branchings. To compare to existing analytical 
results in Fig. [21 we have adopted the same assumption 
here. But the MC algorithm is easily extended to treat 
all incoherent components of a parton shower on equal 
footing and to simulate the physically expected, further 
inelastic interactions of fully formed gluons. In addition. 



we have chosen here the ansatz ^c'^'^^ 



■tqQrg 



CX 



since it is inherent in the existing analytical treatments 
of the LPM effect, and since it serves illustrative purposes 
(we wanted to have in the incoherent limiting case a l/w 
power law, on top of which the characteristic interference- 
induced vti^-modification is identified easily). The pro- 
posed MC algorithm, however, applies equally well to 
arbitrary choices for g-iQT^qQT qt^^]^^ g-qQT^qQTg ^ g^j^^j 

lows us for instance to include a more realistic, steeper uj- 
dependence of a'^Q^^qQrg dose to the kinematic bound- 
ary Ld ~ Eq. Moreover, realistic inelastic cross sections 
provide information about the recoil of the target partons 
Qt- The present proposal may open a new approach to 
follow the dynamical evolution of these recoils. Another 
important challenge of future work is to interface the 
present proposal with the virtuality evolution of a vac- 
uum parton shower. We plan to pursue these questions 
in the context of the MC code JEWEL [13]. 
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